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ABSTRACT 


QOES-VAS satellite retrievals are used to investigate an 
intense oonveotive outbreak over the Mississippi River Valley on 
21-22 July 1982. The primary goals are to assess the strengths 
and weaknesses of three methods for oomputing vertical motion 
using satellite retrievals and to determine the effeots of short 
Interval observations on the calculations. Then, the vertioal 
motions are incorporated with thermodynamic parameters to assess 
the usefulness of VAS data in delineating faotors leading to 
storm formation. 

Results indicate that the quasl-geostrophlc omega equation 
provided patterns and magnitudes most consistent with observed 
weather events and the 12 h radiosonde-derived motions. The vor- 
tioity method generally produced reasonable patterns, especially 
over the oonveotive outbreak, although magnitudes were large due 
to its time derivative. The adiabatic technique gave very poor 
results because diurnal temperature variations extended from the 
surface through 600 mb, and there apparently was a compensating 
effeat in the retrieval algorithm that produoed opposite fluctua- 
tions near 300 mb. These short term variations caused local 
changes to dominant horizontal adveotlon, producing total motions 
that did not agree with synoptic conditions even though patterns 
from adveotion alone were quite good. Thus, results from time 
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derivatives were not as encouraging as hoped, especially for 1 h 
data. Finally, VAS retrievals v/ere of limited use in detecting 
changes in Lifted Indices because fluctuations could be inferred 
from routine surface reports. Although patterns of prooipitable 
water wore reasonable, time variations v/ere suspeot since they 
oould not be confirmed v/ith ground-based data. Nonetheless, in 
general, VAS retrievals were a useful supplement to the tradi- 
tional 12 h radiosonde reports. 
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1 , Introduction 


Synoptio and planotary souloo of motion inoludo many atmoo- 
phorio oyotcmo, but numorouo othor important motoorolegioal 
phonomona ooour on omallor ooaloo. Moaoooalo foaturos havo typi- 
oal wavolengtho of 1-1000 km and poriodo of oovoral houro and 
inoludo oyotcmo ouoh ao thunderstorms, land/ooa br 00200 , and 
frontal ralnbando. Although our knowledge of, and ability to 
forooaot oynoptla ooalo ovonto hao advanced considerably during 
tho past oovoral dooades, corresponding improvements at the 
moooooalo havo boon muoh more limited. Thore aro at loaot two 
explanations for foroaaotlng limitations in mosomotoorology. 
First, mathematical and thoorotical diffioultioo inhibit tho 
development of prediction modols. Second, although largo amounts 
of data are needed to accurately desorlbe these features, only a 
short time is available for data processing. The ourrent 
radiosonde network, having a spacing of 400 Ira and providing data 
at 12 h intervals, is insufficient for this purpose. 

Acourate forecasts of severe conveotivo storms are a major 
challenge for meteorologists. Duo to data limitations, foreoasts 
have been made from a synoptic soale perspective. Soon, however, 
it is hoped that meteorological satellite technology will make a 
major impact on severe storm prediction. Specifically, tho Visi- 
ble Infrared Spin-Soan Radiometer (VISSR) Atmospheric Sounder 
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(VAS), oontalnod on tho moat raoont GOES oorioo, providoo rodi- 
anoo data from which mesoacalo ratriovalo of atmospheric tompora- 
turo and wator vapor oan bo obtainod. Complete doaoriptiono of 
tho VAS inotrumont aro givon by Monzol £t al . (1981), Smith n£ 
nk» (1981b), and Chootoro ofc. aX. (1982). Soundings having opa- 
tial and temporal ooparationo of about 100 km and 1 h aro avail- 
able, mueh amallor than that provided by tho radioeondo not work. 
With thlo now capability, it io hoped that forooaotoro oan uoo 
VAS ooundinge to bottor predlot tho development and ovolution of 
cover 0 thunderstorms. 

Thoro aro throe noooooary ingredients for thundorotorm ini- 
tiation: instability, moisture, and a triggering mechanism* Con- 
cerning stability and moisture content, several recent investiga- 
tions have favorably assessed VAS's ability to detect mososoale 
variations in both time and space (e.g, , Smith at si. , 1981a; Leo 
fife. 1983! Smith, 1983; Chesters ot ni. , 1984). To the 
author’s knowledge, however, there have been no published studios 
to evaluate satellite-derived fields of vortical motion. That is 
tho purpose of tho ourrent paper. 



2. Methods of computing vertical motion 


Vortical motion la an important atmoophorio paramotor 
because it controls the distribution of clouds and precipitation, 
fortioally transports maos, moisture, and momentum, and releases 
the instability of the pro-oonvootivo environment. Since magni- 
tudes of vortioal motion in the free atmosphere are usually only 
several centimotors/sooond for synoptic scalo systems, aocurato 
determination is elusive because diroot physical measurements are 
impossible. Instead, indirect techniques, utilizing other atmos- 
pheric quantities, must bo employed. Miller and Panofsky (1958) 
cite five sohemes for computing large scale vertioal motion, 
i.e., the kincmatio, adiabatic, precipitation, vorticity, and 
numerical (omega equation) methods. All but the precipitation 
procedure are explored in the current study. Eaoh technique is 
based on quite different physical and mathematical assumptions 
and, therefore, possesses different strengths and weaknesses tnat 
are described below. 

K&iejnatte. RQ.tJteA 

The kinematic procedure is based on the continuity equation 
in which divergence of the horizontal wind is integrated upward. 
Vertioal motion in isobaric coordinates ( u> = dp/dt) is given by 

% ° “p+Ap + ^*V 2 >Ap , (1) 
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whore Vg to the horizontal wind and Ap is positive. There are 
□ovoral advantages to the kinematic teohniquo. Thoee inolude 
simplicity of calculation, look of stringent physical assump- 
tions, and the absonoe of time derivatives. On the other hand, a 
significant disadvantage is that orrors in the measured winds 
accumulate during the stop-wise integration to produce unrealis- 
tic at in the upper troposphere. For this reason, adjustment pro- 
cedures, such as those by O'Brien (1970) and Feddar (1981), often 
are used to modify the profiles to a predetermined value (usually 
zero) at the top most level (usually 100 mb). 


Auiabatlo petfrM 


This method is derived from the First Law of Thermodynamics 
by assuming adiabatic flow. Vertical motion can be computed from 


at 


-V 2 «$e - 3 6 / 3 1 
3 0/Tp ’ 


( 2 ) 


where 6 is potential temperature. A simple mathematical formula- 
tion and an insensitivity to wind errors, so that geostrophio 
values can be used, are important advantages. However, three 
disadvantages are a time derivative that may be inadequately 
described by 12 h radiosonde data, inaccurate results when the 
atmospheric lapse rate approaches the adiabatio value, and an 
inability to properly represent dlabatio motions. 
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Si. Ifortloltv mathM 

Quasi-goootrophio thoory io applioablo to largo aoalo sys- 
tems In which thoro Is a quasi-equilibrium between horizontal 
forces. The formulation for the vortioity method is obtained 
from integration of the quasi-geostrophic vortioity equation 
(Holton, 1979) yielding, in isobaric oc rdinates 

“p “ Vap + 17 [ ' H g /3t - V* u e +£),4p • (3) 

where V is goostrophio wind, t is geostrophio vortioity, and f 
B 6 

and f are variable and oonstant values of the Coriolis parame- 
ter, respectively. Strengths of the method are ease of calcula- 
tion, use of geostrophio winds, and the ability to negleot com- 
plex terms such as tilting, friction, and the vertioal advecticn 
of vortioity. Weaknesses inolude possible violation of the 
underlying assumptions during ageostrophlc conditions and, as in 
the kinematic technique, an accumulation of data errors during 
step-wise integration. Finally, as noted for the adiabatic tech- 
nique, the time derivative (of vortioity) may be unrepresentative 
when 12 h data are employed. 

d.. mathal 

The omega equation is derived by eliminating time deriva- 
tives in the vortioity and thermodynamic energy equations. The 
complexity of the result depends on the number of simplifying 
assumptions, ranging from the general ba3 ince equation, 
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containing twelve foroing functions, to the quasi-geostrophic 
omoga equation, having only the two dominant processes of the 
oompioto version. The quaoi-goostrophio omega equation is (Hol- 
ton, 1979) 


(V*+ 




" V 2 [-V g -^TI.<4> 


where I? is the dry air gas oonstant and o (c“51£2), representing 

p09p 

stability, varies only in the vertioal and with time. The left 
side of this equation, a kind of three-dimensional Laplaoian, oan 
be approximated by w multiplied by a negative ooeffieient. The 
right side oontains forcing functions involving the vertioal 
derivative of absolute geostrophic vortioity adveotion and the 
Laplaoian of geostrophic temperature adveotion. Thus, upward 
motion should ooour east of height troughs and in areas of max- 
imum warm air adveotion. A oritioism of this formulation is that 
the two foroing funotions usually are not independent because 
each oontains a common component. Therefore, alternative 

representations have been proposed (e.g., Hoskins si. , 1978; 

Trenborth, 1978). 


Three advantages of the quasi-geostrophic omega equation are 
its explanation of several empirical rules of synoptio meteorol- 
ogy, relative ease of calculation, and a lack of time deriva- 
tives. A disadvantage is the possible violation of quasi- 
geostrophic theory for subsynoptic scale motions and for strongly 
baroolinio systems where diabatic processes are important. 
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Enn£_ /studios 

Based on radiosonde data, the kinouatic method generally 1 b 
considered to provide the most reasonable synoptic scale vortical 
velocities, l.e. , those whioh best correspond to observed map 
features and cloud/preoipitatlon patterns. The numerical method 
appears to be the second ranked procedure. Specifically , Smith 
(1971) found that results from the omega equation were oomparable 
to kinematic values. In a later study, however, the kinematic 
technique w as dearly superior to the omega equation (Vincent efc 
a!* t 1976). This was confirmed by Smith and Lin (1976) who noted 
that the kinematic procedure outperformed both the general bal- 
ance and quasi-geostrophic omega equations. They speculated that 
the kinematic method was better able to capture smaller scale 
processes than the numerical approach. In fact, recent studies 
have shown that mesoscale vertical motions can bo successfully 
calculated using the kinematic technique (Fankhauser, 1969; Moore 
and Fuelberg, 1981; Fuelberg and Printy, 1983). 

Comparisons between different forms of the numerical method 
have produced contrasting conclusions. For example, Smith and 
Lin (1978) noted that the simple omega equation was superior to 
the general balance equation during intense cyclone development. 
However, Krishnamurti and Moxim (1971) and Pagnotti and Bosart 
(1984) emphasized that latent heat release could be an important 
forcing funotion for vertical motion in some situations. More- 
over, Haltiner £t (1963) demonstrated that frictional 
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of foots markedly lnfluonood w fields in the lower troposphere. 

The adiabatio method generally has been a less reliable 
Indicator of largo Boole vertical motion than the kinematic and 
numerical techniques (e.g,, Hansen and Thompson, 1965} Wilson, 
1976). When applied to mesosoale data, adiabatio motions agreed 
favorably with weather patterns in the lower troposphere; how- 
ever, they were inferior to kinematio results in the middle tro- 
posphere (Fuelberg and Lee, 1982). Concerning the quasi- 
gftostrophio vortioity method, Eliasson and Hubert (1953) and Col- 
lins and Kuhn (195*0 observed satisfactory agreement between 
areas of ascent and large scale patterns of precipitation; how- 
ever, its use has been relatively infrequent. 

For radiosonde data, the references cited above suggest that 
each method oan prouuce reasonable estimates of synoptic soale 
vertical motion, with best results from tne kinematic and numeri- 
cal procedures. The current objective is to assess the various 
techniques using VAS satellite retrievals as input. Such an 
evaluation is appropriate beoause satellite and radiosonde data 
are fundamentally different in nature. For example, the 
radiosonde is virtually a point source sensor, while the 
satellite’s radiometer provides volume-averaged measurements, 
thereby achieving smoother definitions of atmospheric structure, 
especially in the vertical. In addition, the satellite’s shorter 
time interval between observations and lack of directly observed 
winds are important differences. Specifically, the current 
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investigation will oaloulate vortioal motions at high temporal 
resolution from the vortioity, adiabatic, and omega equation 
methods and will assess the strengths and weaknesses of each 
technique. The kinematic procedure cannot bo utilized with VA S 
retrievals since only geostrophio or gradient winds can be com- 
puted thermally. The study also will investigate the common 

assumption that a shorter time interval would result in more 
representative vortioal velocities in those methods containing a 
time derivative. Moreover, the ability of each technique to 

bridge the "gap" between 12 h radiosonde data will ue described. 
The final goal is to Incorporate the vertioal motions together 
with satellite-derived fields of water vapor content and stabil- 
ity to diagnose causes for an intense convective outbreak over 
the middle Mississippi River Valley. 



3. Methodology 


a. m aae a daim 

The period 21-22 July 1982 was chosen for study because It 
oontalned significant afternoon conveotion during a period when 
VAS satellite data were available. It is a typiool case in that 
the entire data region wa3 not oloud-free before thunderstorm 
development. Six sounding times were available: 1100, 1300, 
1600, 1700, 2000, and 2300 GMT 21 July. Separate northern and 
southern image sectors of the domain were collected, but, due to 
the small time difference between them (about 30 min), they were 
treated as a single image for retrieval purposes. This consoli- 
dation may be a slight source of error in the 1 h (1600-1700 GMT) 
changes of temperature, vorticity and thermodynamic parameters 
that follow. The 1300 and 1600 GMT image prirs were made during 
the satellited Dwell Image (DI) mode of operation, whiie the 
remaining data were obtained during the Dwell Sounding (DS) mode. 
The DS mode generally is expected to give somewhat better 
retrievals of temperature and moisture because of its larger spin 
budget. Smith .etj, sL, (1981b) give a complete description of 
these modes. 

VAS retrievals were prepared at the NASA/Marshall Space 
Flight Center using techniques incorporating the Man-computer 
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Intoraotivo Data Acooss System (MoIDAS) (Smith fit » 1979). 
Tho physical rotrlovaX algorithm of Smith (1983) was employed; 
tho proooduro addo details of temperature and moisture to a 
"first guess" profile. For tho initial sounding timo (1100 GMT 
21 July), tho first guooo was from a National Meteorological 
Contor ( NMC) 12 h forooaot valid at 1200 GMT 21 July. Aftor 1100 
GMT, VAS retrievals from tho previous timos wore used as first 
guesses. Tho output soundings then wore edited via MoIDAS with 
Inconsistent data being deleted. 

The VAS retrievals underwent a second round of editing at 
Saint Louis University. The objeotives wore to looate any 
remaining anomalous features before objective analysis, obtain a 
working knowledge of the data, and estimate the smallest 
wavelength at which VAS data provided meaningful information. To 
achieve these goals, oonstant pressure analyses of several basio 
and derived parameters were prepared. Results were oheoked 
olosely for vertical stacking, time continuity, appropriateness 
of gradients, and cloud contamination. Deletions were made as 
necessary. Finally, sounding data were created at 100 mb inter- 
vals from the mandatory levels originally supplied. Standard 
temperature interpolation and hypsometric procedures were used. 

ll. Gfr-tecti.ve AnalltBla o£ .MS jcetjcleyala 

Sinoe the convective outbreak occurred over tho central 
one- third of the United States, this was the area of focus. 
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LoQQtlono of edited soundings within the region are shown in Fig. 
1. In selecting tho grid Interval and response characteristics 
for tho objective analysis, tho proooduros of Koch Hi.. 0983) 
wore considered. Table 1 gives tho number of soundings (after 
editing) within tho analysis region at eaoh time. Retrievals 
have an avorago opaoing of 110 km in relatively clear areas, but 
due to oloud oovor over oertaln regions, the separation is not 
uniform (Fig. 1). Therefore, the "random data spacing", the 
average separation for a uniform sounding distribution, was cal- 
culated (Koch Hfc Hi. , 1983). Since the random spacing (Table 1) 
vas always greater than that in dear areas alone (110 1cm), espe- 
cially at 1100 GMT, clustering is indicated. Because of its lim- 
ited and highly clustered data, 1100 GMT was not utilized for 
later calculations. 

Due to the grouping of the data and ambiguity as to the 
smallest scale that VAS soundings are meaningful, meso 0-soale 
resolution was not attempted. Instead, for thermodynamio calcu- 
lations only, the retrievals were treated as a meso a-scale data 
source, assuming an average spacing of 250 1cm. The Barnes ( 1973) 
objective analysis procedure was employed on a grid mesh of 127 
km. Response parameters that retained 62$ of amplitudes at 
wavelengths of 500 km (twice the assumed VAS spacing), i.e., 89$ 
at 800 km (twice the radiosonde spacing), were utilized. Fig. 2 
shows that the objective analysis for geopotential height at 500 
mb for 1600 GMT (middle) has removed most meso 0 -scale features 



Fig. 1. Locations of edited VAS retrievals 
at the six observation times. 
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Table 1. Number of soundings after editing within the 
computational domain. Their separation if 
uniformly distributed, is also shown. 


Sounding 

Time 

(GMT) 

Number of 
soundings 

Random 
data spacing 
(km) 

1100 

81 

245.9 

1300 

127 

191.6 

1600 

155 

171.8 

1700 

120 

197.6 

2000 

157 

170.6 

2300 

170 

163.4 
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Fig. 2. Subjective and objective analyses of VAS-derived heights 
at 500 mb for 1600 GMT 21 July 1982. The top panel is 
the subjective (hand) analysis (SUB) while the middle and 
bottom sections are meso a (MESO) and synoptic (SYNOP) 
scale objective versions. Heights are in decameters. 
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prooont in tho original oubjootivo (hand) analysis (top). Some 
of thooo details nay have boon meteorological phenomena while 
othoro wgpo duo to data uncertainty, but oinoo moot oxhibitod 
poor continuity, thoy v/oro not rotalnod. 

Uhilo mono a-ooalo rooolution wuo appropriate for calculat- 
ing thermodynamic paramotoro, it wao not ouitablo for ppopcring 
vortical motions oinoo oovoral mothodo employed quaoi-gooatrophlo 
thoory. To obtain tho nooooaary larger ccalo fioldo, smoothing 
of tho mooo o-ooalo grido wao oonoidorod, o.g. , Shuman (*957) i 
howovor, it wao decided to begin the analysis proooos anow, this 
time with synoptio scale rooponso parameters that rotalnod only 
32 of amplitudes at 500 km wavelengths (i.o. , 272 at 000 1cm). 
Thus, heights using the new response (bottom of Fig. 2) are con- 
siderably smoother than those at the flnor resolution. Specifi- 
cally, meso a-scalo features are greatly filtered while larger 
scale wavelengths are retained. To verify tho suitability of tho 
analyzed data for quasi-goostrophio calculations, Rossby numbers 
(R Q ) wore oomputed using R 0 a |V 2 |/fL, where L, tho soalo length, 
is one-fourth of tho aotual wavolength. Results of approximately 
0.1 verified the largo scale character of the fields. 

&. fladtfiflaada dafea 

Radiosonde data were obtained from the National Climatic 
Center. Although four times were investigated, i.e., 0000 and 
1200 GMT on 21-22 July 1982, the two middle observations wore 
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moot important oinoo they served as bracketing times for the VAS 
retrievals. Data at 100 mb intervals were prepared using the 
proviouoly mentioned procedures. The Barnes (1973) technique 
again was employed to grid the data, but in this case a 6 0? 
response at 800 km separation (two times the average station 
spacing) was utilized for therraodynamio calculations. As with 
tho VAS data, tho grid interval was 127 km. For oonsistenoy with 
tho satellite-derived vertical motions, corresponding values from 
radiosonde data were calculated from grids again having a 27? 
response at 800 km separation. In addition, due to the availa- 
bility of observed winds, vertical motions also weze calculated 
using the kinematic jrocedure. 



4. Ueathor conditions 


Ao io typical of summer, 21-22 July was not on especially 
baroolinio situation, although there was significant oonvootivo 
aotivity. One tornado, tv;o oases of largo hail, and fifteen 
ooourronoos of damaging Hinds were documented within the grid 
network during the 24 h period beginning at 1200 GMT 21 July. 
Additional severe weather ooourred afterwards. Fig, 3 shows 
National Meteorological Center (NMC) analyses for 0000 GMT 21 
July through 1200 GMT 22 July. It should be noted that the 
poriod contained active shortwave development, not simply advec- 
tlon of pre-existing features. 

At 0000 GMT 21 July, a frontal system extended from the East 
Coast into Kansas. The western half contained only weak gra- 
dients of temperature and dewpoint; however, radar (not shown) 
revealed scattered thunderstorms over eastern Missouri. A cold 
front ourving from North Dakota to Montana was the system of con- 
cern for the middle Mississippi River Valley. At 500 mb, ridging 
extended from Colorado into the Ohio River Valley. To the north, 
a developing shortwave had not yet produced appreciable cyolonio 
flow. No well defined jet streaks were evident at any level (not 
shown) , although strongest winds and the thermal trough were 
upwind of the wave axis, thereby suggesting further oyolonic 
development. 
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3. NMC analyses for the surface and 500 mb between 0000 GMT 
21 July and 1200 GMT 22 July. At the surface, solid 
lines are isobars at 3 mb inter vals while dashed lines 
depict mixing ratios at intervals of 2 g kg“^. At 500 
mb, heights (solid lines) are in decameters while dashed 
lines denote absolute vorticity (10 - ^ s~*). 
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Fig. 3. (Continued) 
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Fig. 3. (Continued) 
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Fig. 3. (Continued) 
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By 1200 GMT 21 July, the weak front from Missouri to 
Colorado had booomo stationary, but the cold front over the Dako- 
tas had progressed southeastward, now extending from Wisconsin to 
Nebraska. Thunderstorm activity was associated with this front 
(not shown); however, the oonveotion over eastern Missouri had 
virtually dissipated. At 500 mb, the shortwave had amplified, 
causing the ridge over the South to become "pinched off" with a 
seoond high over Kentucky. This ridge suppressed convective 
aotivity in the pre-storm environment. A vortioity maximum was 
positioned over central Iowa with positive vortioity adveotion 
(PVA) inferred through Illinois. Maximum winds (not shown) con- 
tinued on the upwind side of the trough, and a jet streak (25-30 
m s -1 ) had developed at 200 mb. 

Renewed storm aotivity over the middle Mississippi River 
Valley, the emphasis of the current study, began near 1735 GMT. 
This was a preferred area for convective development. At the 
surfaoe (Fig. 3), an axis of high mixing ratios ( >16 g kg -1 ) 
extended from the Culf Coast through western Illinois (equivalent 
to dewpoints of about 21 - 23 °C), while at 050 mb (not shown), 
values of 15-16 g kg were evident over western Illinois, In 
addition, temperatures greater than 34°C (not shown) stretched 
from oentral Arkansas into southern Illinois at midday. Thus, 
the region was quite unstable with 1200 GMT Lifted Indiaes (LI) 
ranging from -2 to -5. Given the instability and adequate water 
vapor, only a triggering mechanism was needed to initiate oonvec- 



Hourly satellite imagery beginning at 1600 GOT (Fig. A) 
reveals thunderstorm development in several looations. Although 
largo scale asoont oast of the trough is probable over eaoh area, 
smaller soale forcing appears to be superimposed as well. The 
first area, having tops to 18.9 Ion and extending from 
southeastern Illinois to northern Alabama, likely was associated 
with the nearby gradient of temperature and mixing ratio (Fig. 
3), Although not analyzed by NMC, thiB appears to be a warm 
frontal zone. The second and third storm aroau, with tops to 
16.6 km, were located over Arkansas and southeastern Missouri. 
Based on hourly GOES imagery as well as surface reports, both 
appear to have been induced by separate thunderstorm outflow 
boundaries from earlier convection over western Missouri. In 
addition, surfaoe streamline analyses (not shown) revealed con- 
vergence in these areas. A fourth region of convection developed 
rapidly over southwestern Illinois, just east of St, Louis, as a 
cloud mass progressing across Missouri (Fig. A) moved into a more 
favorable environment that included abundant moisture and large 
scale ascent from the middle tropospheric trough and the advanc- 
ing cold front (Fig. 3). The National Weather Service, recogniz- 
ing these conditions, had issued a severe weather watch for most 
of Illinois just prior to rapid storm development. Tops of the 
ensuing convection exceeded 18.2 Ion. Finally, afternoon heating 
and weak cyclonic flow contributed to thunderstorms along the 
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Fig. 4. Hourly visible satellite iL-igery between 1600 and 2300 GMT 
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Fig. 4. (Continued) 
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Gulf Coast statoo. Of oouroo t tho heating likely onhanoed insta- 
bility in tho othop storm regiono ao well. Fop oxanplo, oatol- 
llto imagery (Fig. 4) shows southwest Illinolo to bo relatively 
oloar before convective development. 

By tho ond of tho otudy period, 0000 GMT 22 July (Fig. 3)i 
tho oold front over tho southeast hod moved llttlo, while tho 
stationary front through Missouri had boon omitted by NMC. Tho 
oold front over tho northwest had oontinued its advanoo, and a 
low had developed along it ovor Illinois in rooponso to tho 
deepening middle- tropospheric system. Thunderstorms wore most 
intense and widespread near this time with strongest oells (Fig. 
4) aligned along the moisture tongue through Illinois and Ken- 
tucky (Fig. 3). NMC now oarriod a squall lino in this region. 
Slight moisture oonvergenoo (reaching 0.7 g kg h ) at 850 mb 
(not shown) was aiding the oonvoctlon. 

At 300 mb, the movement of the deepening shortwave is 
noteworthy. Specifically, it had advanoed into Illinois between 
1200 GMT 21 July and 0000 GMT 22 July, thereby triggering thun- 
derstorm development between the two times. The pattern of 500 
mb absolute vortlcity (Fig. 3) contains a maximum over Wisconsin 
and a lobe extending southward to Alabama. The lobe is ooinoldent 
with maximum convective aotlvity. At 200 mb (not shown), the Jet 
streak remained upwind of the trough. 

Finally, by 1200 GMT 22 July, twelve hours past the primary 
study period, the low had moved to western Indiana, sweeping the 
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oold front through Illinois. HMC now carried a warm front from 
the low toward the oouthoaot. Although, the majority of thunder- 
otoroo had diooipatod, oomo frontal aotivlty oontlnued through 
Indiana. The otrengthoning of tho ourfaoo low wao In rcoponoo to 
oontlnued ohortwavo amplif io jtion at 500 mb (Pig. 3). 



5. Rooulto 


Thio oootion preoonto goopotontial holghto, goootrophio vor- 
tloltloo, ond vortiaol motions ooloulatod from tho radiooondo and 
satellite data. It doooriboo patterns and continuities ot tho 
fivo VAS timos, ao woll ao tholr rolationohipo with tho braokot- 
ing 12 h rndiooondo-dorivod fieldo. For hoighto and vortioitios, 
tho 500 mb lovol will bo examined, but for vertical motlono, 600 
mb is omphaolzod olnoo it mppoared to have tho boot overall 
results. Stability and preoipitablo wator aro tho final parame- 
ters discussed. 

Halfth-t- ani vor-fcloltv 

Beforo doocrlbing the vortical motlono, it is informative to 
present tho more basio parameters of geopotontlal height and 
absolute goootrophio vortioity (Fig. 5)* Continuity among the 
five VAS-dorivod fields is very good at 500 mb in spite of tho 
data gaps noted earlier (Fig. 1). Spoclfioally, the major trough 
through Ivwa and northwest Missouri at 1300 GMT progresses east- 
ward to Wisconsin and Illinois by the last time. The trough 
amplifies through 1700 GMT but appears to diminish thereafter. 
Aotually, however, heights near the trough oontinue to fall 
several meters. The illusion results from strengthening of the 
eastern high pressure ridge (about 50 m) through 1700 GMT, fol- 
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Fig. S. 


Heights (solid, in decameters) and absolute geostrophtc 
vorticities (dashed, in 10"^ s“ l ) at 500 mb. Radiosonde- 
derived analyses are at 1200 GMT 21 July and 0000 GMT 
22 July, while the intervening times are from VAS data. 
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Fig. 5. (Continued) 
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lowed by weakening. The timing of the shortwave in relation to 
the otorm outbroak is noteworthy, /la the trough advances toward 
Illinois and the inhibiting downstream ridge dissipates and 
exits, there apparently is ascent that releases the potential 
instability and initiates the oonvectlvo outbreak near 1735 GMT 
(Pig. 4). This will bo oonflrmod in later seotions. 

Patterns of absolute goostrophio vortioity, computed as the 
Laplaoian of geopotential height, also have very good continuity 
(Fig. 5). Tho main vortioity maximum ( n X n ) associated with the 
shortwave generally propagates slowly northeastward with the 
trough. Values and gradients are greatest at 1700 GMT, reaching 
24 x 10~ 5 s” 1 over western Wisconsin. The elongated minimum from 
the Gulf Coast States through Illinois is well positioned with 
respect to the height ridge, with smallest values ("N") near the 
greatest antioyclcnic curvature. 

A oommon forecasting tool is to qualitatively assess verti- 
cal motion from absolute vortioity adveotion. Based on the omega 
equation (4), upward increasing positive vortioity adveotion 
(PVA) indicates ascent, whereas upward increasing negative adveo- 
tion (NVA) denotes descent. Also, from the vortioity method (3), 
Dutton (1976) and others have noted that when strong PVA occurs 
above approximately 500 mb and/or strong NVA occurs below, there 
is a contribution to negative u. Fig. 5 shows that maximum. PVA 
is associated with the northern portion of the shortwave trough 
and the main cloud area (Fig. 4). Due to orientations of the 
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height and vortioity isopleths, greateet PVA (6.5x10"® s" 2 , not 
shown) and apparent uplift are at 1700 GMT over eastern Iowa and 
southern Wisconsin. Weak PVA enters western Illinois between 
1600-1700 GMT, and is over most of the state by 2000 GMT. This 
adveotion is instrumental in the oonveotivo development. On the 
other hand, NVA, and therefore apparent subsidence, are suggested 
over eastern Kansas and most of Oklahoma, thereby corresponding 
to the relatively clear skies (Fig. 4). 

By comparing patterns from the two data sources (Fig. 5), 
one can evaluate the satellite's ability to bridge the 12 h gap 
between conventional soundings. VAS-derived fields at 1300 GMT 
show good agreement with corresponding radiosonde patterns at 
1200 GMT. The main trough and vortioity oenter move slightly 
eastward during the 1 h period. In addition, the radiosonde data 
verify the pre-storm "capping" ridge and vortioity minimum over 
the middle Mississippi River Valley that is an extension of the 
main ridge over the Great Plains. Similarly close agreements are 
found between satellite-derived fields at 2300 GMT and the 
radiosonde version at 0000 GMT, One again should note the orien- 
tation and slight eastward movements of the trough and associated 
vortioity maximum. Based on the conventional soundings, heights 
in the trough have fallen 15-20 m sinoe 1200 GMT, but as stated 
earlier with the VAS data, this is not obvious from the oontours 
alone because of the downstream ridge. The radiosonde soundings 
also illustrate damping of the ridge over southern Illinois, 
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thereby verifying results from the VAS retrievals. Finally, pro- 
gression of the trough axi* between the two radiosonde times is 
dopioted well by the intervening VAS observations (Fig. 6). For 
eaoh suooessive time, the trough moves further oast by an amount 
proportional to the time interval between observations. The good 
agreements and continuities noted here confirm the operationally 
based findings of Anthony and Hade (1983), 

Although height patterns from the two data souroes agree 
olosely (Fig. 5)» a careful inspection shows that values at 
speolflc locations sometimes exhibit greater discrepancy. This 
will affect results from the vorticity method since it oontains a 
time derivative, but with the numerical procedure only patterns 
at single times are employed. These topics are examined in the 
following sections. 

JtacLLgJJa. method 

Fields of vertical motion from the vorticity method are 
presented for 600 mb in Fig. 7. Radiosonde soundings were util- 
ized at 1200 and 0000 GMT; however, VAS retrievals were input at 
the other times. For both data souroes, backward differenoing 
was employed in the time derivative of geostrophic vortloity 
(35g/3t). An exception was forward differenoing at the first VAS 
time (1300 GMT). Also, for 1700 GMT. instead of utilizing 1 h 
time interval, data at 1300 GMT were employed. Reasons for the 
latter ohoioe are explained later. 



Fig. 6. Approximate locations of trough axes at 500 mb 
from VAS (solid) and radiosonde data (dashed). 
17 denotes 1700 GMT 21 July. 
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Fig. 7. Vertical motions (gb s”*) at 600 mb from the vorticity 
method. Radiosonde soundings were utilized at 1200 and 
0000 GMT while VAS data were employed at the remaining 
times. The dashed line at 2000 GMT is the axis for the 
vertical cross section in Fig. 11. Patterns at 1700 GMT 
are from a 4 h time derivative using 1300 GMT data. 
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Fig. 7. (Continued) 
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Continuity between the five VAS times gone', -ally is very 
good. Areas of greatest upward and downward motion correspond 
oloaoly with the major height and vortioity features {Fig. 5). 
In this comparison* it is important to note that patterns of 
height and vortioity at 500 mb are very similar to those at 600 
mb (not shown). At 1300 and 1600 QMT, the center of ascent along 
the Missouri- Illinois border is positioned ahead of the major 
trough axis. Patterns for these two times appear almost identi- 
cal beoause the choices for time differencing produoe the same 
3?g/3t component* and it is the dominant contributor to u. The 
area of asoent moves little by 1700 GMT, but by 2000 and 2300 
GKT, it propagates into Illinois* coincident with the primary 
storm activity. 

Concerning downward motion* a center extends from southern 
Illinois to western Tennessee through 1700 GMT (Fig. 7), near the 
ridge and vortioity minimum (Fig. 5). By 2000 GMT, however, it 
is replaced by the upward motic'J that initiates the convection. 
This ability of the vortioity method to delineate vertical 
motions over the storm area is very encouraging. A second region 
of positive values is behind the trough over the Northern Plains. 
It shows reasonable continuity by expanding and moving generally 
eastward with the shortwave. A final area of descent, initially 
over Oklahoma and Kansas, is coincident with the ridge and 
minimum vortioity (Fig. 5). Both areas slowly advance northward 
by 2000 GMT, However, by 2300 GMT, the signs reverse, with 
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upward motion now ovor tho Southern Plains. Thin ohongo la duo 
to poaitlvo 35g/3t, ao tho area of minimum vortloity movoo north 
to tho Kanoas-Nebraska border (Fig. 5). 

It la informative to document rolative contributions of tho 
temporal ( (-ac g /at]Ap ) and advootlvo ( f-^g’^Ug + f)]Ap ) com- 
ponents. Fig. 8 shows results at 600 mb Integrated from the sur- 
faoe for 2000 GMT. Tho local derivative (top) is dominant, oven 
though a 3 h Interval is utilized. For individual levels (before 
integration) , the time derivative is most influential below 600 
mb where advection is weak; however, the oomponents are more com- 
parable in the middle and upper troposphere. Nonetheless, after 
integration the upper troposphere (not shown) generally 1 b dom- 
inated by effects of the local derivative from below. It is 
interesting to note that the two components oppose one another in 
the main area of ascent ahead of the trough and in the area of 
subsidence over northeastern Kansas. This is expeoted since 
downwind of troughs, for example, the adveotion term in (3) pro- 
duces desoent in layers containing PVA, whereas the local 
inorease in vortloity produoes ascent. The Inverse is true 
downwind of ridges. Dutton (1976) also noted this tendenoy for 
the terms to oancel one another. In ourrent regions of opposi- 
tion, signs of overall vertioal motion usually are determined by 
the looal derivative. 

Fig. 9 presents vertioal velocities at 1700 GMT based on 1 h 
time differencing from 1600 GMT data. One should recall that the 
diagram in Fig. 7 was based on a 4 h interval. In comparing the 




Fig. 8. Te mporal ( t-H g /3t]flp , top) and advective 

( + f ) ) Ap » bottom) components of 

vertical motion (yb s -1 ) from the vorticity 
method at 600 mb for 2000 GMT . 


40 



Fig. 9. Vertical motions (pb s -1 ) from 
the vorticity method at 600 mb 
for 1700 GMT; 1 h time differ- 
encing from 1600 GMT data was 
employed. 
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two diagrams, sovoral similarities aro apparent, Pirot, both 
voroiono produoo aooont near tho ohortuavo trough and ovor 
woutorn Arkanoao; however, tho primary oontor 1c dioplaood 
further north with tho 1 h difforonolng, oauolng poor continuity 
with tho othor tlmoo. Also, tho 1 h voroion contains upward 
motion ovor wootorn Kanoao that doeo not ooour with tho longer 
time intorval. Hogiono of ouboidonoo from oouthorn Illinois to 
Mississippi, behind the trough ovor South Dakota, and ovor 
Oklahoma and Toxao, are quite olmilar in eaoh oaoe. However, tho 
oontor of descent ovor central Illinois from tho 1 h motiono does 
not appear in the 4 h patterns. Concerning magnitudes, motiono 
from tho 1 h sohemo aro approximately triple the 4 h valuos in 
tho main oontor of ascent over Wisconsin, and three to four times 
greater in tho areas of desoent over Mississippi and Worth 
Dakota* Thus, their agreement with the other times (Pig. 7) is 
poor. At 300 mb (not shown), patterns from both versions aro 
similar, but the 1 h valuos are unroalistioally large, o.g. , 
greatest values aro -68 pb s" 1 and +130 pb s" 1 . 

Differences in the two versions of 1700 QMT vertical motions 
are due to the time derivative and can be explained by the vorti- 
oity patterns of Fig. 5* Specifically, their continuity is quite 
good, but magnitudes at particular locations sometimes vary 
greatly between times. A portion of these fluctuations is due to 
data uncertainty; however, near the trough, suoh errors are 
thought to be a relatively minor component of tho overall change. 



Effooto of duta unoortainty oon bo ovaluatod by introducing 
aimulated orroro into tho obaorvcd goootrophic vortioltloo. To 
maximise tho of foot, motlono duo to tho 3Gg/9t component wore 
oemputod at 1700 GMT from porturbod data at 1600 and 1700 GMT, 
Pig, 10 ohowo rooulto for tho grid point in wootorn Wiooonoin 
that originally had both tho largoot goootrophic vortiolty and 
upward motion. Tho top oootlen olmulatoo vortiolty orroro of +10$ 
(daohod) or -10$ (daohod-dot) of tho oboorvod algobralo dlffor- 
onoo at oaoh lovol. It io obvloua that largo variations in vort- 
ioal motion ooour for a 1 h interval. Of oouroo, with aotual 
data, orroro at Individual lovolo should al tomato signs to somo 
t.ytont, thereby providing partial oanoollatlon upon vortical 
integration. Also, tho offeoto of orroro aro loos if a longer 
time interval io employed. Thio oloarly io indioatod by the bot- 
tom section of Fig. 10 which io baood on tho oamo perturbations 
ao tho top portion, oxoopt there io U h differencing. Those 
rooulto suggest that 1 h time derivatives are not appropriate for 
tho vortiolty method, and that more elaborate procedures may be 
necessary to take advantage of VAS's 1 h capabilities, at least 
for those synoptic situations similar to 21 July 1982. For oases 
with stronger dynomio forcing, the influence of data uncertainty 
should bo smaller; thus, a 1 h time derivative might be more use- 
ful. 

Having established tho continuity of patterns among the five 
VAS times, their agreement with radiosonde-derived fields now is 




Fig. 10. Vortical motions (Mb s”*) from the time derivative com- 
ponent of the vorticity method for 1700 GMT. "O" denotes 
the observed profile for a grid point over western Wis- 
consin. The other two curves represent values based on 
either a +10% (dashed) or -10% (dashed-dot) error in the 
observed algebraic difference in geostrophic vorticities 
at each level. For both panels, 1600 and 1700 GMT data 
were utilised. 
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oonsldored. In goneral, VAS-dorlved pattorno show good con- 
tinuity with those from the traditional aouroo; however » tho 
satellite's values are considerably greater (Fig. 7). At 1200 
GMT, upward motion extends from Minnesota to Nebraska; however, 
at 1300 GMT, tho aentor of aooent over northern Missouri is dis- 
placed considerably. By 0000 GMT, the oentor is along the Xowa- 
Illinols border, and there is a lobe extending southeastward over 
a storm area (Fig. 4). VAS adequately fills tn the 12 h data 
gap, showing the eastward progression of rising motion. Thus, 
oapping desoent over the future storm area at 1200 and 1300 GMT 
ohanges to asoent by 2300 and 0000 GMT. For the subsidence 
behind the trough, there is agreement over the Northern Plains, 
but there are dlsorepanoles in the Southern Plains due to 3e„/3t. 
At 0000 GMT, the descent over Mississippi oorrespor.us to a rela- 
tive laok of oonveotive activity (Fig. 4). A final point con- 
oerns the components of vertloal motion. Although the time 
derivative is dominant for the V AS-derived fields (Fig. 8), the 
terms are comparable in the case of radiosonde-derived motions 
(not shown). One should note that 12 h time differences were 
employed with the RAOB data. Sinoe contributions by geostrophic 
vorticity adveotion are very similar for the two sources, it is 
the greater values of 3? g /3t that lead to enhanced satellite- 
derived vertical motion. 

Vertical stacking of the satellite-derived motions is dep- 
leted in the top oross section of Fig. 11 that passes through the 




Fig. H. Cross section of vertical motion (ub s” 1 ) from the vor- 
ticity method for 2000 GMT. The axis is shown in Fig. 7 
nearby radiosonde sites are indicated. The top profile 
depicts original values while the bottom contains adjus- 
ted values using the O'Brien (1970) technique. 
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convective outbreak at 2000 GMT (see Fig. 7 for axis). Confi- 
gurations at each level generally are simJ'-ar, thereby suggesting 
that horizontal patterns are comparable to those seen at 600 mb 
(Fig. 7). This is verified by results at the individual levels 
(not shown). An unfortunate aspect is that unreailstlo motions 
occur at 100 mb duo to the aooumulation of error upon vertioal 
integration. To eliminate the problem, the 0*Brien (1970) 
adjustment scheme was applied so that values equaled 0 at 100 mb. 
This technique has been used extensively with the lcinematio 
method (e.g., Vincent e£. al. f 1976; Smith and Lin, 1976); how- 
ever, its applicability to the vortieity method has yet to be 
established theoretically. The resulting profile (bottom of Fig, 
11) shows the desired maximum in the middle troposphere. Thus, 
suoh an adjustment scheme, is useful when profiles of vertioal 
motions are desired. Of oourse, it is usually unnecessary near 
the surface and would not be required in the upper troposphere if 
only relative magnitudes are needed at each level. 

& j5aaa^ « g.9.BAcfiB M La .omega eflHa&jaa 

The omega equation (4) includes effects of both temperature 
and vortioity advection. Solutions were obtained using succes- 
sive relaxation with boundary conditions of w c 0 at 1000 mb 
(bottom), 100 mb (top), and on all lateral sides. The constant 
value for the Coriolis parameter was 0.89X10" 11 s" 1 . The lateral 
boundary condition is a potential source of error. Specifically, 
systems moving into or out of the region would produoe non-zero 
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values. In the current case, however, the system of Interest is 
well contained within the computational domain. An additional 
consideration is the stability faotor (a). Values less than zero 
would produoe a hyporbollo formulation with more than one solu- 
tion instead of the usual elliptical result. This problem does 
not ocour in the ourrent study beoause a is always positive. 
Finally, although data errors may produce hyperbollo solutions, 
the large scale response function being utilized eliminates this 
possibility. 

Fig. 12 displays vertical motions at 600 mb from both data 
souj*oes. Continuity between the five VAS-derived fields is 
excellent. Upward motion is consistently ahead of the trough 
(Fig. 5) while there is downward motion behind. The oenter of 
ascent over Iowa at 1300 GMT advances slowly eastward, extending 
from Illinois to northern Arkansas by 2300 GMT. Values are 
greatest at 1700 and 2000 GMT, reaching -2.6 ub s’ 1 over eastern 
Iowa, but magnitudes weaken by 2300 GMT. This variation can be 
explained by the relative contributions of the two equation terms 
[see (it) J . Although they are complementary near 600 mb, PVA 
often dominates the ascent at 1700 and 2000 GMT. Slnae the 
trough is amplifying (Fig. 5), one might expect greatest uplift 
at 2300 GMT. Aotually, however, due to the orientations of 
height and vortiolty, upward increasing PVA is much weaker, 
resulting in the smaller ascent. The oenter of subsidence, due 
mostly bO oold air advectlon, progresses into northwestern Iowa 




Fig. 12. Vertical motions (ub s”*) at 600 mb from the omega equa- 
tion. Radiosonde soundings were utilized at 1200 and 
0000 GMT while VAS data were employed at the remaining 
times. The dashed lines at 2000 GMT are axes for the 
vertical cross sections in Fig. 13. 
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Fig. 12. (Continued) 
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by 2300 GMT, when magnitudes are greatoat. Ovor the Central 
Plains, ascent inoreasoo through 2000 GMT, but the oentor Is vir- 
tually non-existent by 2300 GMT. This feature mostly results 
from warm air advection which diminishes by 2300 GMT; however PVA 
is a major contributor at 2000 GMT. Descent over Oklahoma is 
evident at eaeh time, reaching 1.0 pb s' 1 at 1700 GMT. It 
corresponds to relatively dear skies on the satellite imagery 
(Pig. «). 

Unlike the vortiolty method (Fig. 7), the omega equation 
does not produce downward motion ovor Illinois prior to oonveo- 
tlve development. Instead, weak asoent ooours at each time, 
although values do become more negative as the trough advances 
eastward. Onoe again, it is encouraging that VAS-derlved verti- 
cal motions diagnose the large soale uplift leading to the oon- 
veetive outbreak. 

Major features of VAS-derlved vertioal motions show remark- 
able continuity with those from the 12 h radiosonde data (Fig. 
12). At 1200 GMT, there is weak ascent (-0.5 ub s" 1 ) on the 
border of Iowa, Illinois, and Wisconsin, while weak descent (0.5 
ub a” 1 ) is centered over South Dakota. Satellite-derived pat- 
terns at 1300 GMT show similar positions and magnitudes. By 0000 
GMT, the negative center extends from eastern Illinois to south- 
ern Michigan, while positl e values are over Iowa. Corresponding 
areas at 2300 GMT show close agreement. Finally, weak subsidence 
over Kansas at 0000 GMT supports the sign reversal between the 
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final two VAS times. Although VAS data do a vary good job of 
filling in tho 12 h gap, tho catolli to-derived vortical motions 
generally are largor than those from the aondoo, A possible 
explanation involves the amount of data and its objootive 
analysis. Spooifioally, ovon though response functions were 
identical for both data sources ( 27 % at 800 Ion), tho greater 
number of VAS retrievals may provide bettor resolution of the 
various forcing meohanisms. 

Although there is excellent continuity at 600 mb, it is 
important to consider vertical structure as well. Fig. 13 shows 
two oross sections at 2000 GMT, the first passing through the 
main trough and the second through the storm area (see Fig. 12 
for axes). Vertical stacking is excellent in each case. For the 
northernmost diagram (top), ascent and descent are maximized near 
600 and 300 mb, respectively. In faot, 300 mb always is the 
level of maximum subsidence, with greatest values reaching 3.0 pb 

„i 

s . Thermal advectlon is dominant at this level and generally 
opposes the contribution by vorticlty advectlon. For the south- 
ernmost cross seotion (bottom), regions of uplift associated with 
the trough and with warm air advectlon farther west are clearly 
evident in the lower and middle troposphere. Unlike the vorti- 
olty method, exoessive values do not occur in the upper tropo- 
sphere. The cross seotlons, together with horizontal patterns at 
other altitudes (not shown), indicate that vertical velocities 
are reasonable at all levels. Thut, satellite-derived motions 



Cross sections of vertical motion (jjb s“l) from the 
omega equation for 2000 GMT. Axes are shown in 
Fig. 12; nearby radiosonde sites are indicated. 






from tho quaoi-geootrophio omega equation aro very onoouraging, 
giving fields that have good continuity and that brldgo tho gap 
botwoon radiosonde data. This is espoolally significant in light 
of tho data gaps notod eaplior (Fig. 1). In terms of magnitude, 
they oorroopond to values commonly obtained with NMC forecasting 
models. 

As a final point, vertloal motions wero re-oalnulated using 
a response funotlcn that retained 60? of amplitudes at 800 km 
wavelengths, i.e., tho value employed for radiosonde-derived 
thermodynamic calculations (Sootion 5-g). The goal was to test 
the sensitivity of the vertloal velocities to smaller scale data 
features. Fig. 14 contains resulting patterns of height and 
geostrophio vortloity at 500 mb (top) and vertical motions at 600 
mb (bottom) for 1700 GMT. Compared to the original height fields 
(Fig. 5), the new patterns contain little additional "noise". 
For vortiolcy (Figs. 5. 14), magnitudes are substantially altered 
in the centers, but general patterns and their locations are 
quite similar. The new vortloity field (Fig. 14) only appears to 
contain more horizontal detail than the previous version (Fig. 5) 
because the greater magnitudes lead to an increased number of 
contours. Finally, patterns of vertical motion from the higher 
response (Fig. 14) also are very similar to those in Fig, 12, but 
magnitudes are larger (greatest values at 600 mb are -5.8 yb s" 1 
and 4.1 yb s”^, respectively), and oenters are more evident. A 
review of past studies (e.g. , Barnes, 1984) suggests that these 
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Fig. 14. Heights (decameters) and absolute geostrophic vortlci- 
ties (10 - ^ s _1 ) at 500 mb (top), and numerical vertical 
motions (pb s -1 ) at 600 mb (bottom) for 1700 GMT. The 
input data contained a higher response than that used 
originally in Figs. 5 and 12. 
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lorgor voluoo Bight bo Inappropriate for tho currant oynoptlo 
weather oituatlon. 

Tho Important oonoluoion io that patterns of numorloal vort- 
ical motion are not greatly affooted by tho aholoo of response 
parameters. This finding is similar to that of Stuart (197b) who 
noted that results were not sensitive to tho typo of analysis 
toohnlquo used to propo.ro tho input height data. On tho othor 
hand, patterns from tho vortioity method (not shown) wore greatly 
affooted by tho amount of input detail beoauso of tho time 
derivative and vortloal integration. This relative insensitivity 
to response orltoria is an advantage of tho numerloal prooeduro. 
When ohoosing the parameters, one should oonsidor whether magni- 
tudes of vertioal motion are appropriate. However, it is i por- 
tent to note that higher resolution is more likely to violate 
quasi-goostrophio theory. Therefore, in the current study, the 
response yielding Roseby numbers near 0.1 was utilized. 

1. Adiabatic. jnc-t hod 

Values from this thermodynamically oriented method depend on 
atmospheric stability as well as the local change and horizontal 
adveotion of potential temperature [see (2)]. Vertioal motions 
were obtained between 800 and 200 mb, but not at 900 mb where 
dlabatlo effects usually are quite pronounced. Backward time 
dlfferenolng again was employed to calculate the loos? time 
derivative (36/3t), exoept for forward differenoing at 1300 GMT. 
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Also, olnoo noutral or suporodiabatio lapoo ratoo woro common 
ovor tho wootorn part of the analyolo roGlon botwoon 700 and 500 
mb, potential teraporaturo was required to inoroaoo at loaot 2 
K/100 mb (Fuolborg and Loo, 1982) to avoid unroaliotioally largo 
vertical volooitioo. 


VAS-dorivod vortloal motiona at 600 mb aro givon in Fig. 15. 
Unllko tho prooontation for tho vortioity method (Fig. 7), tho 
1700 GMT map wao obtainod from 1 h difforonoing} a oorrosponding 
dlogrom baood on a A h inoremont will bo ohown lator. Adiabatic 
motions aro not presented for radiosonde data oinoe tho adverse 
VAS-derived patterns to bo dosoribod below would make objootivo 
oomparleons impossible. At 1300 GMT (Fig. 15), there is sub- 
sidence over virtually the entire region. However, asoont then 
begins a general westward progression, covering almost the com- 
plete domain by 230s GMT. In other words, there is a sign rever- 
sal between the initial and final times. These vertioal motions 
aro not associated with tho height patterns (Fig. 5). For exam- 
ple, there is no consistent ascent (desoent) oast (west) of the 
trough and no satisfactory explanation for tho oonveotive out- 
break over tho middle Mississippi River Valley. 


The time derivative causes the unrealistic motions. Com- 
ponents at 1700 GMT (Fig. 16) show that the temporal contribution 
( , top) is at least one order of magnitude larger than 
that from adveotlon ( — , bottom) . Ho meteorologioal pro- 
cess, e.g. , temperature adveotlon, adequately explains the abnor- 



or«m 

OR POOR QUALITY 



Fig. 15. VAS-dorived vertical motions (iib b“*) at 600 mb from 
the adiabatic method. Patterns at 1700 GMT are from 
a I h time derivative using 1600 GMT data. 
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Fig. 16. 


Temporal ( , top) and advective ( , 

bottom) components of vertical motion (yb s - *-) from 
the adiabatic method at 600 mb for 1700 GMT. 
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mally largo middle tropospheric thermal changes oauelng those 
motions. Aftor 1700 GMT, thunderstorms might produoe lou level 
ooollng, and therefore uplift, over the storm region (Fuelberg 
and Printy, 1983), but, of oourse, no oonveotlon ooours in the 
other areas. Thus, problems with 30/3t apparently result from 
diurnal temperature variation?. Normally, diurnal effeots are 
oonflned to the planetary boundary layer, up to about 650 mb. 
However, the satellite retrievals extend them over too deep a 
layer, reaQhing 500 mb at times. These trends are evident in 
Fig* 17 i which is a pressure-time cross section of area averaged 
temperature (T) for three subareas of the total analysis region. 
At 600 mb, there is warming through 2000 GMT in the west (areas 1 
and 2), and through 1600 GMT over the storm region (area 3). For 
example, in subarea 1, T increases 1.23°C at 600 mb between 1600 
and 1700 GMT, and 2.50°C at 800 and 700 mb. For subarea 2, values 
are 1.74°C at 600 mb and Just over 3.00°C at 800 and 700 mb. 
Eaoh area, in fact most of the domain, exhibits ceding 
thereafter. As one expects, the diurnal tendencies are most pro- 
nounoed nearer the surface. 

Based on (2), thermal variability relates directly to verti- 
cal motion. Specifically, increases in temperature (and 9 ) on 
isobaric surfaces result in descent (w>0) from ■ since ae/3p 
is , '.ways negative, while decreases in temperature cause ascent. 
Therefore, over the west, there is subsidence through 2000 GMT 
(Fig. 15), whereas ascent oocars thereafter. Since this is the 




Fig. 17. Pressure-time cross sections of area averaged tempera- 
ture (°C) (following page) for three subareas of the 
total analysis region (this page) . The number in the 
upper right corner of each cross section is the number 
of the subarea on this page. Values are plotted at 
600 and 300 mb. 
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Fig. 17. (Continued) 
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region of largest diurnal fluctuations and most unstable condi- 
tions ( ae/3p) , it is also the area of greatest vertical motions. 
Based on the requirement that 38/3p increase at least 2K/100 mb, 
the 1 h warming of 1 ,7i|°C in subaroa 2 corresponds to average 
desoent of nearly 26 yb . Within the storm region (area 3), 
warming (Fig. 17) leads to subsidence (Fig. 15) through 1600 GMT, 
while cooling produoes asoent thereafter. 

Although the time derivative produoes unrealistic patterns , 
vertical motions frjm the adveotion component alone are very rea- 
sonable at 600 mb (bottom of Fig. 16). (Note that the oontour 
interval is 1 yb a -1 for this map only to Insure pattern recogni- 
tion.) Fields are vary similar to those from the Laplaolan of 
temperature adveotion within the omega equation (not shown). 
Despite these good patterns, the total fields (Fig. 15) generally 
do not correspond to the synoptlo situation (Fig. 5) because of 
exoessive 36/at. 

Based on results from the vortioity method, one might 
presume that adiabatic vertical motions at 1700 GMT from a 4 h 
difference (bottom of Fig. 18) would be superior to those from a 
1 h interval (Fig. 15). Unfortunately, patterns are very similar 
to those at 1600 GMT (Fig. 15) sinoe the dominant time derivative 
was computed fiom 1300 GMT data in each case. Differences 
between the two 1700 GMT maps include opposite signs over Wiscon- 
sin and different locations for the downward motion over the 
west. Neither version is consistent with trough/ridge patterns 




Fig. 18. Vertical motions (yb s~ ) at 600 (bottom) and 300 mb 
(top) for 1700 GMT from the adiabatic method; 4 h 
time differencing from 1300 GMT data was employed. 
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(Fig* 5). Concerning magnitudes, values from the 1 h interval 
generally are more than those from the 4 h version, but this con- 
trast is loss than observed with corresponding fields from the 
vortlolty method (Figs* 7 and 9)* Thus, it appears that neither 
1 h nor 4 h differencing produoos satisfactory results; there- 
fore, oholoe of time intervals is irrelevant. 

The apparently fictitious temperature variations produoe 
corresponding height ohanges (Fig. 5). For example, middle tro- 
pospheric heights near the ridge over the Southeast build in 
response to warming through 1700 GMT, while falls, due to oool- 
ing, occur thereafter. Thus, some of the variability in the 
"oapplng" ridge near the storm area is due to the apparently 
erroneous temperatures; however, the radiosonde data (Fig. 5 ) 
suggest that part of the decay after 1700 GMT is indeed valid. 
Similarly, height ohanges in the western part of the analysis 
region (Fig. 5) are influenced directly by the diurnal tenden- 
cies, However, the shortwave is the primary cause for height 
variations over the Midwest. 

It is important to consider how temperature fluctuations 
affect vertioal motions from the vortlolty method. For horizon- 
tally uniform fictitious thermal tendencies, the heights would be 
uniformly altered. Sinoe 5 is calculated as the Laplaoian of 

o 

height, its patterns and values would be unaffected. Therefore, 
both components of the vorticity method (3) also would be 
unohanged, and the vertical motions would oontinue to represent 
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temperature variations duo to trough/rldgo patterns. On the 
other hand, thermal tondonoles would affoot motions from the adi- 
abatlo method slnoo temperature variability Is considered 
dlreotly In tom 30/3t, In the current oase f It appears that 
most of the temperature fluctuations are either relatively uni- 
form over large portions of the domain (e.g., the West) or are 
dominated by observed meteorologloal phenomena (e.g, * the 
Midwest). Thus, for th* moot part, vortlclty-derived motions 
exhibit reasonable continuity (Fig. 7) and are dearly superior 
to those from the adiabatic method. 

Adiabatic vertloal motions in the upper troposphere differ 
considerably from those at 600 mb. Speoifioally, signs at 300 mb 
(Fig. 19) generally are the reverse of those below (Fig. 15), as 
upward motion at 1300 GMT changes to mostly subsidence by 2300 
GMT. Again notice the laolc of continuity for individual 
features* e.g., the switch from asoent to descent between 1600 
and 1700 GMT over Michigan and Louisiana. The problem, of 
oourse, results from unrealistic 3 0/ 3t. That is, area averaged 
temperatures (Fig. 17) reveal opposite trends between 600 and 300 
mb and the faot that upper tropospherlo thermal variations are 
greatest at this level. For example, early in the period, oool- 
ing ooours at 300 mb over the storm area and two western subre- 
gions, while warming in all three areas and, in fact, the entire 
domain occurs thereafter. Thus, for the reasons stated previ- 
ously, the component of vertical motion due to bs/3t, e.g., 1700 




Fig. 19. VAS-derived vertical motions (yb s - *) at 300 mb from 
the adiabatic method. Patterns at 1700 GMT are from 
a 1 h time derivative using 1600 GMT data. 
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GMT (top of Pig. 20), Is opposite to that observed at 600 mb (top 
of Pig. 16). 

Local temperature ohangos at 300 mb also are not related to 
eynoptlo features (e.g,, horizontal advcotlon). Instead, thoy 
apparently result from tho physical retrieval algorithm. Specif- 
ically, this procedure appears to compensate for excessively deep 
diurnal ef feats by produoing tho reverse situation in tho upper 
levels (V. L. Smith, personal oommunioation, 1985). In the few 
areas where low level diurnal variations are small and extend 
over only a shallow layer, there is little or no upper level com- 
pensation. Although upper tropospherio fluctuations obviously 
affeot goopotentiid heights at 300 mb (not shown), low level 
thermal variations usually dominate. Nonetheless, vortioity- 
derlved motions at 300 mb (not shown) again are greatly superior 
to the adiabatlo versions. Contributions by thermal advection 
(bottom of Fig. 20) still generally are smaller than those from 
3 0/3t (top), but advection is more important than at 600 mb (Fig. 
16), especially at 2000 QMT. Compared to the thermal advection 
term of the omega equation (not shown), patterns are quite simi- 
lar; however, values from -V 2 * V0 are larger. 

Finally, 300 mb vertical motions at 1700 QMT based on a 4 h 
time interval (top of Fig. 18) are no better than the 1 h ver- 
sion, For example, the 1 h field (Fig. 19) oontains large sub- 
sidence over Miohigan due to warming; however, the 4 h map shows 
asoent over the same region. Thus, thermal trends over one time 



ORIGSMA* ' 

OF 'POOH Qu.' «. 



Fig. 20. Temporal 


( ’ top) and advective < = mt!% 

bottom) components of vertical motion (pb s“^) from 
the adiabatic method at 300 mb for 1700 GMT. 
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interval oan bo quite different from thooo for a longor period, 
and unfortunately, the adiabatlo teohniquo ie very aonaitivo to 
thooo variations. Neither the 1 h nor 4 h versions is Qocoptubio 
□inoo patterns do not oorrospond to hoight configurations (Pig. 
S). An attempt was cade to ralnimizo the "diurnal effect" by oal- 
oulating 1700 GMT vertical motions from data at 1300 and 2300 GMT 
(10 h interval). Unfortunately, the time derivative still is 
dominant. Thus, results at 600 and 300 mb (not shown) are only 
sxightly improved. 

In summary, because of the time derivative, the adiabatic 
method does not appear to be a promising teohniquo for computing 
vertical motion from VAS retriovals. Patterns exhibit little 
relation to synoptlo features because implied diabetic effects 
violate the underlying assumptions of the method. Patterns from 
the advcotion component are favorable, and that forolng function 
is incorporated into the omega equation whose results were quite 
aooeptable. 

&. Klncmfrlg. jnafthpA 

Past studies of vortical motion (see Section 2-e) generally 
have considered the kinematio method to be a "standard" against 
which other procedures are oempared. Unfortunately, however, it 
cannot be applied to satellite-derived data since only geos- 
trophic or gradient winds oan be derived thermally. Neverthe- 
less, it is informative to briefly describe general similarities 



botwoon radiooondc-dorivcd kinematic motions and versions that 
are avail ablo from VA8. 

Pig. 21 ohowo kinonatic motions at 600 mb. At 1200 GMT 21 
July, thoro lo doooont in the pro-convective anvironmont ovor 
southern Illlnolo, and uplift, aaoooiatod with tho trough (Pig. 
6) , ovor eaetorn Nebraska, By 0000 GMT 22 July, tho suboldonoo 
moves oaotward ao accent advances into Illinoio and Indiana. This 
area extends Into Alabama, thoroby onoompaoeing tho ontlro region 
of oonvootlvo activity (Pig. A). Also, an aroa of doooont 
devolopo ovor Missouri behind tho trough. Idontioal reoponso 
parameters (27? at 800 km) were used to analyze both tho 
radiosonde and satellite data; however the kinematic technique is 
not limited by quasi-geostrophio theory. Therefore, vortical 
motions also wore calculated using a higher response (602 at 800 
km). Resulting magnitudes (not shown) are somewhat larger than 
those in Pig. 21 , but patterns are very Bimiler. 

Radiosonde-derived motions from the vorticity, numerical, 
and klnematlo techniques should be oompared briefly (Figs. 7, 12, 
21). At 1200 GMT, the vortioity (Fig. 7) and numerical (Fig. 12) 
procedures produoe areas of asoendlng and descending motion asso- 
ciated with the trough, while the klnematlo method (Fig. 21) 
shows only the area of asoent. The centers from eaoh technique 
are displaced from those of the others. Over the future storm 
region, there is subsidence from all but the omega equation. At 
0000 GMT, all three schemes show good agreement for the upward 
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Fig, 21. Radiosonde-derived vertical motions (yb s" 1 ) 
at 600 mb from the kinematic method for 1200 
GMT 21 July and 0000 GMT 22 July. 
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motion ahood of tho ohortwavo (Pigs. 7» 12, 21); however, asoont 
from tho vortioity method is displaced farther west. In addi- 
tion, oaoh deplots desoont behind tho trough, although in 
slightly different locations, and asoont over tho cjnveation. 
Concerning magnitudes, results from tho vortioity and Iclnematio 
techniques are comparable but are larger than those from tho 
omega bqnatlon. In summary, oaoh version is related to the main 
height featuros (Fig. 5)f moreover, no one prooeduro is dearly 
superior, 

sctalls about the individual techniques have been presented 
in earlier sections; however, it is informative to briefly com- 
pare basic features of the Iclnematio and VAS-derived vertical 
motions. Only the 1300 and 2300 GMT satellite observations are 
considered since they are closest to the radiosonde times (1200 
and 0000 GMT) , and adiabatic motions are not inoluded since 
results already were shown to be quite poor. With the vortioity 
method, patterns (Fig. 7) show general agreement with those from 
the kinematio technique (Fig. 21). Specifically, at the initial 
times, both versions deplot subsidence near southern Illinois and 
ascent near the trough which, by the final times, also extends 
over the storm activity (Fig. 4). However, at the end of the 
period, the downward motion behind the trough shows less agree- 
ment, and values are much larger with the vortioity method. 

Vertical motions from the omega equation also exhibit simi- 
larities with those from the kinematic technique. Specifically, 
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the main areas of ascent (Figs. 12 and 21), whioh show only fair 
agreement at the Initial times, have muoh oloser correspondence 
later In the day. In addition, at the end of the period, both 
contain aesoent behind the trough, although in somewhat different 
locations. Finally, a dissimilarity at 1200-1300 GMT is the area 
of descending motion that the omega equation does not depict over 
southern Illinois. 

VAS and klnematio vertical motions were oompared beoause 
kinematic values generally are very reliable. Nonetheless, even 
that procedure has limitations (see Seotion 2-a). Therefore, the 
satellite-based technique that agrees olosest with klnematio 
results is not necessarily superior. Instead, the optimum method 
will yield vertioal motions that 1 ) correspond most with the 
height field, 2) show the greatest horizontal and vertioal con- 
tinuity, and 3) best relate and fill in the 12 h gap between 
radiosonde-derived motions using the same prooedure. 

£. fignBaclaga sL jggfeftgfla 

Previous sections qualitatively assessed the individual 
methods for obtaining vertical motion. This seotion compares the 
techniques to one another using both qualitative and quantitative 
procedures. The first evaluation is to average vertical motions 
in three key regions: 1) behind, and 2) ahead of the main trough, 
and 3) over the primary storm area (Fig. 22) . "Boxes" 1 and 2 
propagated eastward with the trough; however, "box" 3 remained 




Fig. 22. Areas for which average vertical motions were calcula- 
ted. At the top, box 1 is behind the shortwave trough, 
while at the bottom, box 2 is in advance of the wave. 
Box 3 over southern Illinois is centered over the storm 
outbreak. Boxes 1 and 2 propagate with the trough al- 
though box 3 remains fixed. 13 denotes 1300 G11T. 
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stationary. The Lagrangian approach for the northernmost boxes 
depiots the evolution of vertical motions near the developing 
shortwave. The Eulerian approaoh for the storm area shows 
motions prior to and during the oonveotive outbreak. Pressure- 
time cross seotions for the vortloity, numerical, and adiabatic 
methods are presented for boxes 1-3. The reader is encouraged to 
review the previous space maps at 600 mb to aid in the comparison 
(l.e», Pigs. 7, 12, 15). 

Behind the trough (Box 1, Pig. 23), the vortloity and 
quasi-geostrophio omega equation provide the most reasonable pat- 
terns. (Note that the contour interval is 0.5 lib s" 1 for the 
omega equation bat 2.0 pb s” 1 for the other methods). Specifi- 
cally, the vortloity method depicts a gradual transition from 
ascent to descent during the 10 h period of shortwave amplifica- 
tion, with greatest subsidence at 2300 GMT corresponding to max- 
imum development. Of oourse, the cross aeotion would be more 
appealing if an adjustment scheme, suoh as O'Brien (1970), had 
been utilized as in Pig. 11. Por the omega equation, there is 
descent in the middle and upper troposphere throughout the 
period, with greatest values at the final two times. The adia- 
batic technique offers the least satisfactory patterns. The 
"diurnal and compensating" effeots disoussed previously are 
dearly evident, and there is no trend toward increasing downward 
motion behind the shortwave. 
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East of the trough (Box 2, Fig. 24), the omega equation pro- 
vides the best results sinoe it consistently shows upward motion 
that is maximized in the lower and middle troposphere at 2000 
GMT. The vortlolty method also provides good continuity, except 
at 2000 GMT. Although horizontal depletions reveal some ascent 
east of the trough (Fig. 7)» the placement of the box at this 
time mostly oaptures the larger area of subsidence. An encourag- 
ing aspect of the oross section is that values generally peak 
between 500 and 300 mb without using an adjustment technique. 
Patterns from the adiabatic technique again are dominated by tem- 
poral changes in potential temperature. Therefore, the low level 
asoent after 1700 GMT is not attributable to the trough. 

Finally, over the storm area (Box 3, Fig. 25), the vortlolty 
method yields sneouraging patterns. There is descent through 
1700 GMT with an abrupt change to asoent thereafter. This 
result is consistent with the height fields (Fig. 5), the 
radiosonde-derived vortlolty and kinematic motions (Figs. 7, 21), 
and the timing of the convective outbreak (Fig. 4). Moreover, 
greatest uplift ooours at 2300 GMT, near the time of maximum oon- 
veotive intensity (Fig. 4). The numerical method also performs 
well. Although there is no subsidence prior to storm develop- 
ment, weak asoent is maximized in the lower levels at the final 
two times. Therefore, like the vorticity method, the trend is 
toward increasing upward motion. Adiabatic values vividly show 
the effects of diurnal ohanges in the lower levels and resulting 




of the shortwave trough) . 
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compensation aloft. Reoalllng the horizontal maps (Figs. 15. 
19) » the oross seotlonal patterns do not relate to the oonveotlve 
development t despite that appearance In this depletion. As an 
alternate representation, Fig. 26 contains vertioal profiles from 
the time oross sections In Fig. 25. Notioe the abrupt ahange in 
vertioal velocity from the vortloity method, the gradual 
Inoreases with the numerical teohnique, and the variations In the 
adiabatlo scheme attributable to dlabatio prooesses. 

In addition to qualitative comparisons, two quantitative 
assessments of spaoe-tirae continuity are provided (Chesters 
.al. . 1984): 1) the RMS soale length for gradients within the 
field at time t, and 2) the temporal correlation between succes- 
sive fields. 

The gradient soale length, L(t), is given by 

L ( t ) * [<AV 2 (t)>/ <[dV(t)/dx - <dV(t)/dx>] 2 >] ?2 , (5) 

where 7 represents vertioal motion at a grid point, AV(t) is 
given by AV(t) a V(t) - <V(t)>, and ensemble averages <...> are 
taken horizontally over the entire field. Large L(t) indicates 
smooth patterns, whereas reduced values depict smaller scale 
features. 

Scale lengths for gradients of vertical motion at 600 and 
300 mb are given in Fig. 27. Also inoluded in the left panel are 
reBultB for absolute geostrophio vortloity at 500 mb (Fig. 5) 
and, in both panels, the contribution of potential temperature 
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advoofcion In tho adiabatic toohnlque. At 600 mb, patterns from 
tho quasi-goootrophio omoga oquatlon generally aro cmoothor than 
thoao from tho vortloity method. Thio aontraot lo moat pro- 
nounood at tho two tlaoo utilizing radiosonde data (1200 and 0000 
GMT) when tho dlfforonoo avoragoo 116 tern. Tho amount of dotall 
provided by olthor aoherao dooo not ohango greatly during thfc 12 h 
period, a foot readily verified by tho horizontal dopiatlono 
(Pigs. 7 » IS). An Interesting point la that scale lengths of 
VAS-derlved goootrophlo vortloity are somewhat greater (smoother 
patterns) than results from the omoga equation or vortloity 
method. This is expected slnoe derivatives of vortloity influ- 
ence the vertical motions. Finally* adiabatic motions are by for 
the smoothest. Soale lengths of its advcotlon component are more 
similar to those of the other parameters, thereby confirming that 
total fields are controlled by the time derivative. 

Results at 300 mb (Fig. 27) generally are similar to those 
at 600 mb. For example, patterns from the omega equation again 
generally are smoother than those from the vortloity method. A 
marked difference, however, is that soale lengths for the adia- 
batlo method are considerably smaller than at 600 mb (except at 
2300 GMT), whereas they are larger for temperature advection. 
The greater detail In upper level adiabatic motions, especially 
at 1700 and 2000 GMT, results from tho Increasing Importance of 
horizontal adveotlon. Nonetheless, earlier sections demonstrated 
that patterns from this procedure were not related to synoptic 
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foaturoo. 

Ad a cooond of fort to quantitatively aososo continuity, tem- 
poral correlation ooofflalonto wore computed. Tho defining equa- 
tion io (Choatero skill,, 1904) 

C(t.t') ° <AV(t)AVCt')>/[<AV 2 (t)><AV 2 (t , )>] Ia . (6) 

whore t and t* are tho two observation times. In general , high 
correlations lndloate a continuity of patterns; however, results 
oan bo misleading. For example, meteorologically significant but 
rapidly propagating features aan produce small or oven negative 
correlations, whereas slower propagation in an unrealistic direc- 
tion oan yield large values. To further complicate matters, the 
time intervals being tested are not constant but range from t to 
3 h. Thus, tho coefficients have limited usefulness, but unfor- 
tunately, no other statistic appears to provide a more meaningful 
evaluation. 

Figure 28 presents correlation coefficients for 600 and 300 
mb. Oeostrophio vortloity at 500 mb (bottom panel) shows con- 
sistently high correlations, indicating good temporal continuity. 
Thus, a major input to the vortloity and numerical schemes shows 
a reasonable propagation of features (also see Fig. 5). With the 
vortloity method, two of the values at 600 mb are misleading. 
For the 1300-1600 GMT comparison, the value 0.99 is due to the 
domlnanoe of the time derivative component, whioh, as noted ear- 
lier, was the same in eaah oase. On the other hand, the negative 
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Fig. 28. Temporal correlation coefficients of vertical motion at 
600 and 300 mb and absolute geostrophic vorticity ("Z") 
at 500 mb (bottom). Abbreviations are as in Fig, 27. 
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value for 1700-2000 GMT rofloota rapidly propagating but oontinu- 
oua horizontal pattorno (Fig. 7) • Tho omoga equation generally 
oxhibito high oorrolationo at 600 mb, except for 2000-2300 0MT 
(oee Section 5**o for reaoona). For both the vortioity and numer- 
ical methods, there are low correlations between 1200 GMT 
radiosonde- and 1300 GMT VAS-derlved patterns, but values are 
considerably higher at the 2300-0000 GMT pairing, when the 
shortwave is better developed and more fully oontolned within the 
computational region. Correlations for the adiabatic method ini- 
tially are high at * 00 mb, but this is misleading since they 
represent good continuity among erroneous features. The value 
-0.88 between 2000-2300 GMT, ref loots the sign reversal in verti- 
cal motion. Results for temperature adveotion range between 0.24 
and 0.85. Correlations at 300 mb (Fig. 28) generally are similar 
to those below; however, a notable exception is the 2000-2300 GMT 
value from the omega equation. Moreover, adiabatio motions, 
although generally unrealistic, have a higher correlation between 
2000-2300 GMT because the adveotion component is as influential 
as the time derivative. 

&. lh9rji<?.dy. n PHl<?. fa ctor-?. 

As stated in Seotion 1 , three ingredients are required for 
thunderstorm formation: instability, moisture, and a triggering 
meohanism. For the oonveotive outbreak of 21 July, the shortwave 
trough has been shown to produoe large scale asoent over the mid- 
dle Mississippi River Valley. The temporal and spatial evolution 
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of tho two thermodynamic factors now lo considered. 

Flguro 29 presents radiooondo-dorived Lifted Indioes (LI) at 
1200 GMT 21 July and 0000 GMT 22 July, as well as V AS-derived 
Indioes from 1600-23 JO GMT, Information concerning the response 
oharaoteristios of the input data was given in Sections 3-h and 
o. For the special times, instead of using surface data provided 
with the VAS retrievals, high density hourly reports from the 
airways network were employed. Thus, the calculations utilized 
the greatest amount of surfaoe data possible. Indioes at 1200 
and 0000 GMT were oomputed according to the standard definitions; 
however, values from 1600-2300 GMT were based on a teohnique 
developed at the National Severe Storms Forecast Center (Hales 
and Doswell, 1982). In this proorJure, mean boundary layer data 
are not used; instead, surface parameters are assumed to 
represent the entire layer during that part of the day when it is 
well mixed. Since this assumption would not be valid In the 
early morning, "surface LIs" were not calculated at 1300 GMT. 

At 1200 GMT (Fig. 29), strongly stable conditions are over 
Michigan while lower values cover the midsection of the United 
States, especially along the Missourl-Illinois border. By 1600 
GMT, the entire analysis region apparently is less stable; how- 
ever, most patterns are similar to those seen earlier. An excep- 
tion is the lobe of relative stability over Missouri which is due 
to a pocket of drier air at the surface (not shown). Values con- 
tinue to decrease over the Mississippi River Valley through 2000 
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Fig. 29. Radiosonde- (1200 and 0000 GMT) and VAS-derived Lifted 
Indices. For all times, surface data were from ground 
based observations. 
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GMT, roQohing a minimum loss than -10 over southwost Illinois and 
northern Mississippi. Growing Instability also is depleted over 
the Central Plains. By 2300 GMT' oondltlons begin to stabilize 
from Illinois to Mississippi' although Indices oontinue to 
deorease over Kansas and Nebraska where oonveotion begins after 
the study period. Finally' at 0000 GMT, conditions appear to be 
more stable than 1 h earlier. When Interpreting VAS-derived sta- 
bility indices, it is Important to emphasize patterns and trends 
rather than specific' values. This is neoessary because a very 
unstable index based on radiosonde data may be a common 
ocourrenoe with VAS retrievals, thereby suggesting that a dif- 
ferent set of criteria must be employed. Further researoh is 
needed to establish these values. Anthony and Wade (1983) 
observed that VAS-derlved stability parameters were less stable 
than sonde-derived versions although patterns and trends were 
oomparable. 

The trend toward increasingly negative Lis over the Missis- 
sippi River Valley through 2000 GMT certainly supports the oon- 
veotive outbreak. However, an interesting question is the pro- 
portion of these variations that is due to the hourly surface 
data versus that part attributable to the satellite's 500 mb tem- 
peratures. To answer this question, new Lifted Indices were com- 
puted from 1600-2300 GMT (Fig. 30). Although the surface data 
were identical to those employed previously, radiosonde-derived 
500 mb temperatures for 1200 GMT were used at each of the four 



Fig. 30. Lifted Indices using 1200 GMT radiosonde-derived 500 

mb temperatures at each time along with hourly surface 
reports. 
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times, Therefore, differences in trends between the two versions 
of LI will be attributable only to the changing 500 mb tempera- 
tures provided by VAS. 

Patterns of LI in Fig. 30 are very similar to those in Fig. 
29. For example, the most stable air is over the northeast, 
relatively unstable conditions are over the lower and middle Mis- 
sissippi River Valley and Central Plains, and a pooket of rela- 
tively stable air is in oentral Missouri. An important point at 
eaoh time is that values in Fig. 29 are approximately 2 units 
less than those from the new version (Fig. 30). This difference 
is oaused by VAS's 500 mb temperatures being consistently colder 
than the 1200 GMT sonde data. 

To examine the stability fluctuations, changes in LI (ALI) 
were computed between successive times for each set (Fig. 31). 
Patterns are very similar, e.g. , each depicts negative ALI over 
the storm area through 2000 GMT and increases thereafter. Also, 
the expansion and enhancement of instability over Kansas and 
Nebraska through 2300 GMT is found in both versions. Thus, the 
Important conclusion is that VAS's 500 mb data provided little 
additional information about variations In LI during this oase. 
That is, the sonde's unchanging 500 mb temperatures together with 
hourly surfaoe reports were sufficient to describe the evolving 
patterns. It must be emphasized that VAS retrievals would be 
valuable in those situations having rapid variations in middle 
tropospheric temperature. Finally, other stability parameters 
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Fig. 31. Changes in Lifted Indices with varying VAS-derived 500 mb 
temperatures (top) from Fig. 29 aud constant radiosonde- 
derived 500 mb temperatures (bottom) from Fig, 30. 17-20 

denotes 1700 to 2000 GMT. 
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wore computed for this oase as well, o.g. , K Indox and Total 
Totals; however, alnoe their patterns were similar to those of 
LI, they will not be dlsoussed hero. Their agreement indicates 
that trends at 850 mb were similar to those at the aurfaoe. Of 
course, VAS data would be useful in oaloulating these lndloes 
when conditions at 650 mb were different from surface reports 
and/or underwent fluctuations between the 12 h radiosonde obser- 
vations. 

The third requirement for thunderstorm initiation is suffi- 
cient moisture. Horizontal analyses of surface-300 mb preoipit- 
able water from sonde data at 1200 and 0000 GMT and from VAS 
retrievals between 1300-2300 GMT are presented in Fig, 32, At 
1200 GMT, a moist region extends from Minnesota to the Gulf 
Coast, and VAS-derived patterns at 1300 and 1600 GMT generally 
continue this feature. Between 1600 and 1700 GMT, just prior to 
the oonveotive outbreak (Fig. A), moisture content increases 
greatly (10 mm) over the lower Mississippi River Valley (Fig. 
33). However, from 1700 to 2000 GMT (Figs. 32-33) , there is a 
corresponding decrease over the some region that continues 
through 2300 GMT (Fig. 32). As noted for stability, patterns and 
trends of preoipitable water probably should be emphasized rather 
than speoiflc values. 

Although the values depict a favorable trend in re! allot to 
the storm outbreak, tin variations are quite large. Hillger 
( 198^1) noted h changes in precipltable water reaching 11.6 mm, 
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Fig. 32. Radiosonde- (1200 and 0000 GMT) and VAS-derived pre- 
cipitable water (mm) for the surface to 300 mb layer. 
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Fig. 32. (Continued) 


96 


ORIGINAL. PAC : i;i 
OK POOR QUALITY 



Fig. 33. Changes In VAS-derived precipitable 
water (mm) for the surface to 300 mb 
layer, 16-17 denotes 1600 to 1700 GMT. 
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and ho argued thoir validity on tha baoio of low lovol moiaturo 
adveotion and ourfaoo moiaturo oonvorgonoo. Smith nt £l. 
(1981a) oboorvod 3 h fluotuationo noar 30 mm. Zn tho currant 
otudy, thoro la insufficient moiaturo advootlon or moiaturo oon- 
vorgonco to oxplaln tho 1 h lnorcaaoa of 10 mm. Moroovor, our- 
faoo dowpointa, about (not shown), remain noarly oonotant 
from 1600-1700 QMT. Similarly, moiaturo tranoport and ourfaoo 
dowpolnt fluotuationo do not oubotantlato tho dooroaooo in pro- 
olpitablo wator between 1700-2000 OUT (Pig. 33) • Thoroforo, tho 
variations must bo oonoldored ouapoot, and they may result from 
tho rotrloval algorithm. Recalling Seotion 5-d, warming below 
600 mb (and subsequent height rioao) ooourred over the aroo. 
through 1700 GMT, with ooollng (and height falls) thereafter. 
Those ohangoa were not oonflned bo tho boundary layer. Further- 
more, as a result of the deep, unrealistic warming, the retrieval 
algorithm may have produced inoreaoing preolpltable wator beforo 
tho oonveotlon with tho reverse occurring thereafter (W. L. 
Smith, personal aommunloatlon, 1985). However, over the west, in 
spite of the diurnal temperature fluctuations, there is little 
variation in preolpltable water although that area was much drier 
originally. Additional research is needed to clarify this rela- 
tionship. Nonetheless, oven though the abrupt changes appear 
ouapoot, the overall moist tongue extending along tho Mississippi 
River Valley is quite real since it is verified by radiosonde- 
derived fields and surfaoo data. 
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Tho oquoo for tho oonvootlvo outbroal: now lo oloar. Tho 
ohortwovo trough oauaod upward motion thQt roloaocd tho onorgy of 
tho preferentially moiot and unotablo air ovor tho Miooioolppi 
Dlvor Valloy loading to thundorotorm dovolopmont noar 1735 QMT 
(Fig* 4). VAS data rovoal tho throe roquiremonto in thio area, 
although thoro aro imporfeotiono in accurately diagnoaing oach 
paramotor. 



6. Summary and conclusions 


A oaso atudy hao investigated tho intense oonvootivo 
activity that ooourred over the middle Mississippi River Valley 
on 21-22 July 1982. Tho goal woo to aoaooo tho relative 
otrengtha and weaknesses of varioua prooedurea for calculating 
vertical motion from VAS aatollito retrievals. Magnitudoa and 
oontinuitieo of motiona from the 1-3 h soundings were evaluated 
along with their ability to fill in the 12 h gap between 
radiosonde data. In addition, the study considered the 
satellite's performance in diagnosing those variations in stabil- 
ity and humidity that were conducive storm formation. The 
research benefitted from having VAS retrievals at 1300, 1600, 
1700, 2000, and 2300 GMT 21 July. This allowed an examination of 
conditions prior to and during the thunderstorm outbreak that 
began shortly after 1700 GMT. For thermodynamic calculations, 
VAS retrievals were objectively analyzed to achieve meso a -scale 
resolution of atmospheric structures; however, for vertloal 
motloits, only synoptic scale resolution was permitted* 

Using the satellite-derived input, the quasi-geostrophic 
omega equation generally produoed the best overall results. 
These motions exhibited excellent spatial and temporal continuity 
and consistency with the observed weather. Moreover, they prop- 
erly bridged the 12 h Interval between radiosonde-derived 
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motions. Tho vortioity method provided the oeoond beet results. 
Although patterns noar tho etorm area wore reasonable, their gen- 
eral continuity and rolation to weather events and sonde-derivod 
motions were less fovorablo. The adlabatio teohniquo produoed by 
far tho worst results, Patterns showed poor continuity and had 
little or no oorrespondenoe with synoptic features. Finally, 
sonde-derived kinematic motions, whloh served as a general com- 
parison for VAS fields, agreed reasonably with weather condi- 
tions. 

Poor results from the adiabatio teohniquo occurred because 
diurnal temperature variations through 600 mb violated the adia- 
batic assumption. Furthermore, tho retrieval algorithm appeared 
to compensate for the exaggerated diabatio tendencies by produc- 
ing the opposite effects in the upper troposphere. Thus, verti- 
cal motions showed unrealistic patterns and magnitudes because of 
the time derivative term. Excessively large magnitudes were 
influenced by the short 1-3 h interval between observations, and 
values would have been even greater if unrealistic neutral or 
ouporadiabatio lapse rates between 700-500 mb had not been 
suppressed. Adlabatio vertical motions depended directly on 
thermal fluctuations on individual isobario surfaoes. Therefore, 
lower tropospheric motions from the time derivative often had no 
oorrespondenoe with those in the upper levels. On the other 
hand, since the adveotion component was not based on temperature 
changes between times, it produced very reasonable patterns. 
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Results from tho vortiolty method wore muoh bettor than adi- 
abatlo values beoauso vertloal motions depended on the Laplaolan 
of geopotontial height that was oaloulated from vertically 
Integrated temperatures. Therefore, uniform height variations 
would not affeot the vertloal motions. Tn the ourront case, 
local vortlolty ohanges wore dominated by che presence of the 
trough/ridgo patterns. Since tho contribution by adveotion also 
was linked to the shortwave, overall vertloal motions were 
appropriate, although magnitudes were quite large, due to the 
small interval between observations. 

Tho omega equation contains no time derivative; Instead 
vortical motions depend only on fields of temperature and geos- 
trophic vorticity at single times. Since these patterns ware 
very reasonable, both components of u showed good horizontal con- 
tinuity. Therefore, the shortcomings of the other two techniques 
did not apply, and the best overall results were obtained. 

Before computations began, it was hoped that the short time 
interval between satellite observations would improve vertical 
motions in those methods containing a time derivative. Unfor- 
tunately, however, results from the time components were disap- 
pointing. Temporal variations in potential temperature did not 
appear to be meteorologically significant, i.e., they seemed 
mostly attributable to limitations in the data. Changes in vor- 
ticity were more reasonable. Thus, for cases such as 21 July, 
having data at 1 h intervals apparently will only be useful in 
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tho adiabatic and vortioity teohnlquoo when data uncertainties 
are roduood. In situations of strongop foroing, l.o., a stronger 
meteorological signal, inoreaood oloud cover likely would limit 
tho number and quality of soundings. Since the numerloal method 
is not as greatly af footed by data errors, it appears to be the 
preferred teohnique for oaloulating vertioal motions at this 
time. Its use allows one to take best advantage of the short 
interval soundings provided by VAS without results being overly 
influenced by data ambiguities. 

Concerning the thermodynamic parameters, the increased fre- 
quency of upper air information was not valuable in diagnosing 
variations in the Lifted Index sinoe fluctuations on 21 July were 
dominated by surface effects that oould be obtained from conven- 
tional data. Finally, patterns of VAS-derived preolpitable water 
were similar to those from radiosonde data. The VAS data did 
reveal rapid ohanges in magnitudes, but these could not be con- 
firmed on the basis of surface reports. 

During this case, VAS provided useful information to supple- 
ment the usual 12 h radiosonde soundings. However, only a single 
ease was evaluated, and to the author's knowledge, it was the 
first to evaluate VAS-derived vertioal motions. Therefore, many 
other periods having a variety of synoptio conditions need to be 
examined. In addition, a new physical retrieval algorithm has 
been developed in which excessive diurnal and compensating ‘f,em- 
perature variations are reduced (W. L. Smith, personal communica- 
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tion, 1965). Thus, an ospooially Intones ting study would bo to 
repeat the ourront oaloulations based on tho Improved soundings. 
Finally, bettor procedures for handling data gaps should be 
developed. In the ourront oaso, despite numerous gaps, rosults 
did not appear to bo greatly affected. Only through suoh future 
studies can VilS's ability to provide tho mososoale data needed 
for improved forecasting be clearly evaluated. 
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